Filling the disk hollow following binary black hole merger: 
The transient accretion afterglow 
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Tidal torques from a binary black hole (BHBH) empty out the central regions in any circumbinary 
gaseous accretion disk. The balance between tidal torques and viscosity maintain the inner edge 
of the disk at a radius r ~ 1.5a — 2a, where a is the binary semimajor axis. Eventually, the 
inspiraling binary decouples from disk and merges, leaving behind a central hollow ("donut hole") 
in the disk orbiting the remnant black hole. We present a simple, time-dependent, Newtonian 
calculation that follows the secular (viscous) evolution of the disk as it fills up the hollow down to 
the black hole innermost stable circular orbit and then relaxes to stationary equilibrium. We use 
our model to calculate the electromagnetic radiation ("afterglow") spectrum emitted during this 
transient accretion epoch. Observing the temporal increase in the total electromagnetic flux and 
the hardening of the spectrum as the donut hole fills may help confirm a BHBH merger detected by 
a gravitational wave interferometer. We show how the very existence of the initial hollow can lead 
to super-Eddington accretion during this secular phase if the rate is not very far below Eddington 
prior to decoupling. Our model, though highly idealized, may be useful in establishing some of the 
key parameters, thermal emission features and scalings that characterize this transient, ft can serve 
as a guide in the design and calibration of future radiation-magnetohydrodynamic simulations in 
general relativity. 



I. INTRODUCTION 



Binary black holes (BHBHs) are among the most 
promising sources of gravitational waves detectable 
by gravitation- wave detectors such as LIGO [l|, 
VIRGO 0, il, GEO [|, and TAMA [1, 0, as well as 
by the proposed space-based interferometers LISA 
BBO ^] and DECIGO [13]. The development of stable 
algorithms to integrate Einstein's field equations of gen- 
eral relativity numeric ally in 3-1-1 dimensions, such as the 
BSSN formalism (Til [T2| and the generalized harmonic 
approach [l^, together with the identification of suit- 
able gauge conditions, enabled several pioneering simu- 
lations that demonstrated how to track the late inspiral, 
plunge and and merger of a binary black hole (BHBH) 
in vacuum [ij, [iM EB] • More refined follow-up simula- 
tions of these strong-field, late phases, joined onto to 
analytic, post-Newtonian calculations of the early in- 
spiral epoch p^ . are now cable of producing accurate 
gravitational waveforms for merging BHBH binaries with 
companions spanning a range of mass ratios and spins. 
These theoretical waveforms will be important as tem- 
plates both for the detection of BHBH binaries and for 
the determination of their physical parameters, such as 
the masses and spins of the binary companions and the 
black hole remnant. 

With the binary BHBH merger problem in vacuum 
well in hand, it is now important to turn to the problem 
of BHBH coalescence in astrophysically realistic gaseous 
enviromcnts in general relativity. When the orbit of the 
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BHBH binary is sufficiently close, the ambient gas should 
have little influence on the BHBH dynamics, but the im- 
pact of the black hole on the gas will be considerable. 
In particular, the capture, shock heating and accretion 
of gas may result in appreciable electromagnetic radia- 
tion. Such an "afterglow" could provide another observ- 
able signature of BHBH coalescence in addition to gravi- 
tational waves. Electromagnetic radiation can also serve 
as a useful probe of the gas in galaxy cores or in other 
environments where mergers may take place, as well as a 
probe of the physics of black hole accretion. 

The effect of a BHBH system on a rotationally sup- 
ported, geometrically thin, circumbinary disk has already 
triggered several important investigations, most in New- 
tonian gravitation. One key finding is that the Bardeen 
& Petterson mechanism 18| will align the disk to the or- 
bital plane of the binary on a viscous timescale [l^, . 
Another result, first appreciated in the context of disks 
around pre-main-sequence binaries and protoplanets or- 
biting a central star, is that the disk is truncated at an 
inner edge rodgc where gravitational tidal torques aris- 
ing from the binary and viscous stresses in the disk bal- 
ance each other [2l|. As a result, a hoUowed-out region 
("donut hole") forms in the disk, causing gas accretion 
onto the black holes to remain well below the rate it 
would have in the absence of tidal torques. During the 
binary inspiral, viscosity in the disk keeps the ratio of 
the inner edge radius to the semimajor axis a of the bi- 
nary roughly constant: r edge /a ^ 1.5 — 2 [2ll. [2^. [23l. [23| . 
Eventually, the binary separation shrinks sufficiently that 
the binary inspiral timescale due to gravitational wave 
emission becomes shorter than the viscous timescale in 
the disk. At this point, the binary "decouples" from the 
disk, and the inspiral proceeds all the way to p lunge and 
merger while the disk barely evolves [l^, Ha 113 ■ The 
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black hole remnant quickly settles down to a quasista- 
tionary state, after which the disk gas begins to drift 
inward on a secular, viscous timescale, filling the hollow 
down to the innermost stable circular orbit (ISCO) of the 
remnant black hole at risco- Eventually the disk settles 
down into quasistationary equilibrium and accretes at a 
steady rate onto the remnant. The scenario described 
here is summarized schematically in Fig. [TJ 

In this paper we follow the secular viscous evolution of 
a gaseous accretion disk following decoupling. We per- 
form a simple, Newtonian calculation to determine the 
time-dependent inspiral and structure of a geometrically 
thin disk as it fills up the hollow region and flows into 
the black hole remnant. We also calculate the associated 
thermal radiation spectrum of the disk from the time 
of decoupling to the the establishment of disk equilib- 
rium. Our model is highly idealized in that we assume 
a radiation spectrum at each radius given by thermal 
black body emission at the local effective temperature. 
In a numerical example chosen to illustrate the model we 
also adopt a constant (turbulent) viscosity in space and 
time. While idealized and largely pedagogical, the model 
provides a reasonable first approximation to the relevant 
time and length scales characterizing such a scenario, as 
well as the typical magnitudes and frequencies of the as- 
sociated accretion-driven electromagnetic afterglow. The 
model may be useful in designing and calibrating future, 
more detailed, simulations that employ more realistic mi- 
crophysics and relativistic gravitation. Finally, adopting 
such a simple description for this transient yields a calcu- 
lation that exhibits considerable scale freedom. In partic- 
ular, the time-dependent solution to the nondimensional 
evolution equation that governs the gas inflow is uniquely 
determined by specifying only two free parameters in the 
initial data, the radius of the disk at decoupling, ri and 
the BH remnant ISCO, risco- All other parameters may 
be scaled out of the problem. Moreover, given these two 
parameters, the equation does not depend on any further 
details of the disk microphysics. 

Several recent calculations and hydrodynamic simu- 
lations have studied the prompt response of the cir- 
cumbinary gas disk to gravitational wave-driven BHBH 
mass loss and remnant recoil following merger (see, e.g., 
[H m [13, iH, m and references therein). The simu- 
lations treat effects that can occur on a hydrodynami- 
cal timescale in the absence of kinematic (shear) viscos- 
ity. They arise from those perturbations of the disk that 
are induced by the time- varying gravitational field of the 
binary just prior to merger and the recoiling remnant 
after merger. In fact, viscosity is typically ignored in 
these simulations. By contrast, the scenario considered 
below focusses on the slower, secular release of energy 
that ultimately accompanies the viscosity-driven inflow 
and accretion of the bound, rotating gas in the disk onto 
the remnant. The idealized model we adopt strictly ap- 
plies only when the recoil speed of the central remnant 
is low, in which case the remant remains near the cen- 
ter of the hollowed disk as it accretes. For nonspinning 



black holes, the recoil speed is identically zero for equal- 
mass companions and small for companions of very un- 
equal mass, falling roughly as q^, where 5 < 1 is the 
mass ratio. It reaches a maximum of about 175 km/s 
for q « 0.36 (see [33| for a review and references). For 
spinning black holes of comparable mass the recoil can 
be much larger, but it is very sensitive to the spin mag- 
nitudes and orientations |34|] and again decreases as q^. 
Recoil speeds are estimated to be low (< 200 km/s) in 
most galactic mergers [sHi . Recent simulations of BHBHs 
orbiting inside circumnuclear disks in merged galaxy rem- 
nants give a median recoil speed that is quite low, < 70 
km/s 36]. Determining the distribution of pre-merger 
BHBH masses, mass ratios, spins and spin orientations 
remains an area of active investigation. We anticipate 
that in generic cases, both hydrodynamical (e.g shock) 
and secular (viscous) dissipation of energy are generated 
in the disk following decoupling and contribute to the 
radiation afterglow. We also anticipate that turbulent 
magnetic fields will provide the main source of (effective) 
viscosity in a realistic accretion disk. Future simulations 
will therefore need to incorporate magnetic fields to treat 
all phases of of the afterglow properly. 

We adopt geometrized units and set G = 1 = c below. 



II. BASIC MODEL 

We focus on the epoch following the decoupling of 
the binary BHBH from the circumbinary accretion disk. 
Our calculation begins when the inner radius of the disk 
reaches the decoupling radius ri and we follow its evolu- 
tion as it fills in the hollow around the merged remnant, 
M. We adopt Newtonian physics throughout. We as- 
sume that just prior to decoupling the balance between 
tidal torques from the binary and viscous stresses in the 
disk results in a near-equilibrium disk which has the pro- 
file of a stationary thin disk with an "effective ISCO" at 
ri. We take the outer boundary of the disk to be infinite. 
This constitutes our initial data. 

As the gas loses angular momentum and diffuses in- 
ward from ri following decoupling, we assume that it 
continues to inspiral in Keplerian circular orbits about 
the central black hole in the thin-disk approximation. 
We adopt axisymmetry about the rotation axis of the 
disk and neglect the self-gravity of gas in the disk. We 
take the inner boundary of the inspiraling disk to be the 
remnant black hole ISCO, risco, where we assume that 
viscous torques vanish. 

In our model we assume that the energy dissipated by 
turbulent viscosity in the disk is all emitted as thermal 
black body radiation at the local effective temperature. 
Thus the net flux is a time-dependent superposition of 
local black body spectra from concentric annuli in the 
evolving disk, both inside and outside ri. The magni- 
tude of the integrated (total) radiated flux does not de- 
pend on the black body assumption and is independent of 
the detailed microphysics governing the (vertical) struc- 
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FIG. 1: Schematic figure of a Keplerian accretion disk orbiting a coalescing binary black hole, (a) The balance between tidal 
torques and viscosity maintains the inner edge of the disk at redge ~ 1.5a — 2a, where a is the binary semimajor axis, (b) 
Eventually the binary reaches a critical separation at which its inspiral timescale due to gravitational wave emission equals 
the viscous timescale in the disk, (c) The disk then decouples from the binary, which coalesces, leaving a central black hole 
remnant. Viscosity in the disk drives gas toward the remnant, filling up the hollow between the decoupling radius ri and the 
ISCO radius risco. 



ture of the disk, including its local scale height, viscosity, 
pressure, temperature, and opacity. In an illustrative ex- 
ample we set the turbulent viscosity v to be constant, 
noting that in some a-disk models the viscosity is found 
to be a weak function of radius [23|. The magnitude of 
the viscosity, which does depend on the adopted micro- 
physics, establishes the physical timescale for the inward 
diffusion of gas in the disk, but it can be scaled out of 
the numerical evolution of the disk and the emission. The 
disk microphysics is also important for determining stan- 
dard "modifications" to the black body spectrum arising 
from, e.g., the dominance of electron scattering over true 
absorption in the inner region of the disk, as well as from 
comptonization or synchrotron radiation [37| . These ef- 
fects tend to increase the color temperature above the ef- 
fective temperature of the radiation and thereby harden 
the emitted spectrum. We postpone any consideration of 
these refinements to a future analysis. 



III. BASIC EQUATIONS 

A. Disk Structure 

The evolution equation for the surface density E(i,r) 
in an axisymmetric, geometrically-thin disk may be de- 
rived by combining the equation of mass conservation 
(continuity equation) and the equation of angular mo- 
mentum conservation (the (^-component of the Navier- 
Stokes equation). The result is [33| 



dt r dr \ dr \ 



(1) 



where r is the cylindrical radius measured from the disk 
axis (or black hole remnant), t is the time and v is the 



viscosity, which can be quite general at this point (i.e., 
V — i>{t,r)). Equation ([T|) assumes that the angular ve- 
locity of the disk is Keplerian, Q. = (M/r^)^/^. Taking 
the accretion rate at infinity to be independent of time, 
we seek to solve the equation subject to the following 
boundary conditions: 



b.c.'s 



0, 



constant, r oo 



(2) 



For initial conditions, we shall take the system im- 
mediately after black hole merger at the center of the 
disk hollow. At merger the disk has not changed appre- 
ciably from its structure at decoupling, given that the 
viscous timescale is much longer than the BHBH inspi- 
ral timescale after decoupling. Prior to decoupling the 
balance between tidal torques and viscous stresses pro- 
duces a near-equilibrium disk with an "effective ISCO" 
at ri — Todgo l-5a — 2a, where a is the semimajor axis 
of the binary. Consequently, we take the initial disk to 
be an equilibrium solution to equation ([T]) (i.e. a solution 
to dYi/dt — 0) satisfying the boundary conditions 



i = : i/E = 



(i/S])oo, r ^ oo 
0, r = ri 



(3) 



The initial profile of the disk is therefore a "donut" de- 
scribed by 



i = : vT.= 



^ { (i.S)oo (l - {ri/rf'^) , r>n 



r < ri 



(4) 



In fact, the profile of the final equilibrium disk, formed 
once the "donut" hole fills with gas down to risco and 
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settles down, has a similar, familiar form, 



B. Afterglow Radiation 



t = oo 



0, 



r < 



(5) 

The goal is to solve equation llj for the transient evolu- 
tion during the interval < i < oo. 

The accretion rate M{t,r) — 2TTrJ^Vr, where Vr is the 
inward radial velocity of the gas, may be calculated from 



M{t, r) 



dr 

i/2 5(rV2^E 



dG_ 

dr ' 



dr 



(6) 



In equation ([6]) the quantity G{t,r) is the viscous 
torque fssj . 



G{t,r) = 2nr^iyY. 



dr ' 



(7) 



evaluated for Keplerian motion. Using equation ^ to 
evaluate the accretion rate as r oo yields 



Mr. 



Af(i,oo) = 3tt{vJ:)c 



constant 



(8) 



Equation ^ relates the asymptotic accretion rate to the 
asymptotic disk structure, both of which remain constant 
during the transient. This scenario is consistent with the 
notion that conditions far away from the central binary 
control the rate at which gas is fed into the disk. 

To solve equation ([1]), which in general is nonlinear, it 
is convenient to introduce the following nondimensional 
variables: 



(r/ri)i/2, S^S/Soo, y = st, 

v/^oo, T = t/Uis , (9) 



where 



Ari 
3 i^c^ 



(10) 



is a viscous timescale. In terms of these variables, equa- 
tions ([T]) becomes 



dy 1 d^ {yv) 
d^ " 



(11) 



^2 Qg2 ' 

and must be integrated subject to boundary conditions 

s, s ^ oo 



b.c.'s : y = 
and initial conditions 

t = 0: y = 



0, s = (nsco/?'i)^^^ 



(S-1)/P, 5> 1 

0, s < 1 



(12) 



(13) 



The rate of viscous dissipation per unit surface area of 
the disk is given by [lH 



, , 9 M 



(14) 



Evolution equation ([T]) or (|lip must be solved to deter- 
mine the transient dissipation rate (jl4p as the disk fills up 
the hollow and then causes the region outside the original 
hollow to adjust in response. Assuming that this energy 
is all radiated locally, and approximating the emission 
to be a thermal blackbody radiation, the local disk sur- 
face temperature Ts{t, r) may be equated to the effective 
temperature, which is determined by the dissipation rate 
according to 

aTt{t,r) ^ D{t,r) ^^lyY.^ , < t < oo , (15) 
8 r-^ 

where a is the Stefan-Boltzmann constant. For the initial 
and final equilibrium disks described by equations @ and 
([5]) the result can be found immediately: 



where 



(3MMc 



nr^) 1 - (roq/r) 



1/2 



0, 



' cq 



^cq ^isco 



r > r, 

r < To, 
(16) 



(17) 



oq 



For the transient, the temperature must be determined 
numerically by integrating equation (jlip and substituting 
the result into equation p3|) . 

Given the surface temperature, the transient specific 
flux F^{t) measured by an observer at distance d whose 
line of sight makes an angle i to the normal to the disk 
plane is determined by integrating over the entire disk 
surface, 



2tt cos i 

d2 



B,{Ts{t',r))rdr, (18) 



where B^{Ts{t' ^r)) is the Planck function, t' = t — d is 
retarded time and v is the photon frequency [soj . Equa- 
tion (jl8p is best evaluated in terms of a nondimensional 
function f*{t' , x) of nondimensional frequency x accord- 
ing to 



_ 27r cos n S aT^ ^ , 

fvy-)— ^2 _5 ,, ''isco/ V- -i^j-i 



where 



r(i',^) 



du U ; 

1 exp(a;T*/Ts 



1 



(19) 



(20) 



and where we have introduced the parameters 



MlMo^/^TTrl 



5 



(21) 



The quantity provides an estimate of the characteris- 
tic temperature in the main radiating region of the final 
equilibrium disk, which resides near rjsco, and hv,, is the 
characteristic frequency of the emitted thermal radiation 
from this region. The specific luminosity L^it), summing 
over both surfaces of the disk, is related to Fi, [t) accord- 
ing to 



L,{t) = -F,{t). 

cos I 



(22) 



The total luminosity L{£) integrated over all frequencies 
is then given by 

L{t)= dvL^it) = -aTtrl^l r{t' ,x)dx {23) 
Jo ~ Jo 

The the initial and final equilibrium disks provide ana- 
lytic temperature profiles (i.e. equation I16p . from which 
the fluxes can be determined without knowing the tran- 
sient evolution. At low frequencies x ^ (nsco/foq)'^^^, 
where r^q is defined in equation (I17p . we have 



r(teq,:E)«:rr(8/3)C(8/3)xi/3, 



(24) 



where r(8/3) = 1.50457 and C(8/3) = 1.28419. To derive 
this result we have used the approximation T^,/Ts{tcc[,r) ~ 
{r/r-iscaf'^^ in the outer disk where the low- frequency ra- 
diation is generated, together with the identity 



r(s)C(s) 



exp(a::) — 1 



dx. 



(25) 



Equations pQ]) and ([24]) combine to give the well-known 
low-frequency component of the thin-disk thermal spec- 
trum, Fi, '-^ v^l^ . The total integrated luminosity and 
fiux from the initial and final disk also can be obtained 
analytically. Using the identity 



f*{t' x)dx 



45 r, 



cq 



m equation (|23p gives 



4 _4 2 ^isco 



(26) 



(27) 



or, equivalently. 



Equation ([28)1 recovers the standard result for the total 
luminosity of a thin accretion disk with an inner bound- 
ary at Toq. This result is customarily obtained directly 
from the integral 

L{t) = 2 / aTf 2TTrdr = 2 / g(i^S)— 2nrdr, (29) 



substituting equations ([4]) and ([5]) for the two equilibrium 
profiles. Its significance for afterglow radiation is that, 
at least in our simple model, the total electromagnetic 
fiux will increase from its value immediately after BHBH 
merger (as signaled by peak gravitational wave emission) 
to a value that is larger by a factor of ri /risco as the disk 
settles into final equilibrium. The timescale over which 
the increase will occur will be the viscous timescale near 
the decoupling radius in the disk, ri . In the next section 
we calculate the time-varying evolution and electromag- 
netic emission during this transition. 



IV. THE TRANSIENT AFTERGLOW: AN 
ILLUSTRATIVE EXAMPLE 

A. Nondimensional Description 

Equation (|11[) is a parabolic, nonlinear, diffusion equa- 
tion with a variable diffusion coeffient (e.g., the coeffi- 
cient decreases as for the linear case with constant 
D). In general, such a nonlinear equation must be inte- 
grated numerically, which, in anticipation of future ap- 
plications, motivates our approach here. However, there 
are standard analytic (e.g. Green function) techniques 
that can be implemented for the linear case We 
adopt a Crank-Nicholson finite-difference scheme, which 
is fully implicit and second order in both space and time. 
Such an implicit scheme is well-suited to accommodate 
the large spatial extent of the disk that must be covered 
to implement the boundary conditions. This coverage is 
facilitated by introducing a logarithmic grid in s from 
s = (nsco/fi)^/^ < 1 to s > 1 dO]. Since a Crank- 
Nicholson scheme is unconditionally stable for numerical 
timesteps of arbitrary size, the step size is restricted only 
by accuracy considerations and not grid spacing and can 
be much larger than that required by the usual Courant 
condition for explicit schemes. 

To illustrate the transient behavior with a simple 
model, we take the viscosity in the disk to be a constant, 
in which case iy — I. Equation (jlip is then uniquely 
determined by setting rjsco and ri. We set rjsco = 6M, 
the value appropriate for a Schwarzschild black hole rem- 
nant [4l|, and ri = lOOM, a choice consistent with some 
recent estimates of disk properties at binary BHBH-disk 
decoupling [42]. 

The growth of the surface density E(t, r) inside and 
outside the hollow following decoupling and merger is 
shown in Fig. [2l The associated growth in the local ac- 
cretion rate M{t,r) (see equation [S]) is plotted in Fig. [31 
The secular transition from the initial to final equilib- 
rium disk is evident as gas migrates inward toward the 
black hole remnant. 

The transient thermal spectrum (equations [T9l and [20l) 
is plotted in Fig. [H The spectrum hardens and the total 
flux and luminosity increase in magnitude as the hollow 
flUs with (radiating) matter. In analogy to equation (j2ip . 
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FIG. 2: Profiles of the disk surface density S at se- 
lect times T during accretion following the binary BHBH 
merger. The disk decoupling radius is at ri = lOOM 
and the ISCO is at risco = 6M, where M is the mass 
of the black hole remnat. The initial (final) equilibrium 
disk is indicated by the heavy solid (dashed) curve. Pro- 
files at intermediate times are shown by dotted curves at 
r = 0.0694, 0.208, 0.347, 0.694, 1.39, 2.78, 6.94, 69.4 and 6940. 
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FIG. 3: Profiles of the local disk accretion rate M at select 
times T following the binary BHBH merger (see equation [6}. 
The curves are labeled as in Fig. [2] 



FIG. 4: The frequency distribution of the emitted electromag- 
netic thermal radiation at select times r following the binary 
BHBH merger. The nondimensional distribution function /* 
and photon frequency x are defined in equations (|19p - (|21|l . 
The curves are labeled as in Fig. [2] 



define the quantities 



hvi — kTi 



(30) 



which characterize conditions near the inner edge of the 
initial disk. The peak frequency evolves from near vi 
to over the course of the transient, where v^/iyi — 
(ri/risco)'^^* ~ 8.2 in our example. Over the same time 
frame, the total luminosity increases according to equa- 
tion whereby L{oo)/L{0) = ri/nsco ~ 16.7 in this 
example [43| . 

The increase in total luminosity is plotted as a func- 
tion of time in Fig. [5] The characteristic timescale over 
which the transient lasts is St ^ 10 — 100, or, accord- 
ing to equations ® and Uni), (5t - (10 - 100) ivis- We 
note that the timescale ivis defined by equation ^ is 
related to the viscous timescale at the inner edge of ini- 
tial disk, ri, where the initial surface density, l](0,ri), 
is identically zero (equation 0]) . The gas that establishes 
the final equilibrium flow at the ISCO originates from 
radii r >> ri, where the the density S(0,r) is nonzero 
and closer to its final equilibrium value. Since the viscous 
timescale on which this gas drifts inward varies like r^, 
the equilibration timescale is therefore much longer than 

^vis ■ 

The luminosity increase with time is directly correlated 
with the increase in the accretion rate at the ISCO, as is 
evident from Fig. [51 That the two curves plotted in the 
figure line up so closely once the transient gets underway 
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FIG. 5: The increase in luminosity L and ISCO accretion rate 
Misco with time r during accretion following the binary BHBH 
merger. Note that the time coordinate for Misco represents 
retarded time. 



this end we define several auxiliary quantities, includ- 
ing the black hole mass parameter Afg — M/10^Mq, 
the distance parameter dio = d/10 Mpc, and the ac- 
cretion rate in Eddington units, ttie = M /Me, where 
Me = 12Le/c^ = 1.68 x lO^^Mg gm/s is the Edding- 
ton accretion rate for the final equilibrium disk [i^l and 
Le = 1.26 X 10'*'* Me erg/s is the Eddington luminosity. 
In terms of these parameters, we have 



" ^ * ' ISCO 

L{oo) 



0.906 X 10^ mJ^'^M"^''^ K, 
1.89 X 10*6 m]/V/g"*/^ Hz, 



■70 1 • 1/4 7i,r-l/4 

78.1 A'lg ' 



eV/h, 



1.67 X 10-24 rn^^Ml'^d-^ 



erg 



1.26 X 10''* TTiE erg/s. 



*" cm2 s Hz' 

(34) 



The above quantities set the scale for the characteristic 
innermost disk temperature, peak thermal photon fre- 
quency, observed thermal flux and total thermal lumi- 
nosity for the final equilibrium disk with risco = 6Af. 
The fourth line provides the scale factor that appears in 
equation ([TO]) for the specific flux. Using this factor to- 
gether with Fig. |4] yields the the transient spectrum in 
physical units. 

The corresponding parameters for the initial disk at 
decoupling, with ri = lOOAf , are given by 



(r 0.1) can be explained as follows: during this epoch, 
the surface density in the main radiating, inner region is 
found numerically to satisfy the approximate relation 

v^it, r) « C{t) (l - (n,eo A)*/') , (31) 

where the factor C{t) climbs from near zero initially to its 
late, asymptotic value C(oo) = [vYj)^ (see equation [SJ . 
Using equation ([5T|) in equation (^5)) yields 



L{t)/L{oo) « C{t')/{vJ:), 



(32) 



Using equation in equation ^ to calculate Afisco(i)j 
and inserting equation ([5]), then gives 



Afi3co(i')M^oc«i(i)/i(oo), 



(33) 



which explains the numerical result revealed in Fig. O 
Consistent with the fit (|3ip for the surface density is the 
implication that the accretion rate in innermost region 
near the ISCO must be nearly independent of r during 
the transient, M{t,r) « 37rC(t). This result is verified 
numerically in Fig. [31 



B. Physical Units and Scaling 

To be astrophysically useful and applicable to general 
cases we need to restore physical units to our nondimen- 
sional numerical solution and provide appropriate scal- 
ing with the independent physical parameters. Toward 



m 



0.653 X 10^ V/g K, 



1.36 X 10*^ 



m'J^M-'/' Hz, 



5.63 m^/^Afg-*/'' eV/h, 

0.695 X 10-24 m^^'Hll'^d^^ , 
^ ^ *" cm2 s Hz 

0.754 X 10*^ niE erg/s. 



(35) 



The fact that the thermal emission hardens and increases 
in magnitude during the transient is again evident by 
comparing equations ([M)) and ([55)1 . For black holes of 
mass Afg « 1 accreting at itie ~ 0.1 the peak emission 
at the start of the transient lies in the UV band and 
then hardens to EUV and soft X-ray radiation at late 
times 45]. 

None of the physical disk structure or emission fre- 
quencies and magnitudes calculated above depend on any 
details of the disk microphysics in our idealized model. 
However, we do need to probe the microphysics in order 
to assign a physical timescale to the transient. The di- 
mensional time parameter r can be translated into phys- 
ical time only when ivis appearing in equation (jlOp has 
been evaluated, and for that we need to know the viscos- 
ity. Our numerical implementation assumed a constant 
viscosity in the disk. To get a realistic estimate of the ap- 
propriate value to use we can adopt the standard Shakura 
and Sunyaev [4^ a-disk model and employ the value it 
gives for the turbulent viscosity at ri. In this model, 
the radius ri will reside in the "inner region" , where 
the pressure is radiation-dominated and the opacity is 
electron scattering-dominated, provided it lies within a 
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critical radius rjn , where [13, 13 

njM « 2.56 X 10^ a2/2i^W2i^2/2i_ ^gg^ 

For those cases in which ri < rin [23| the viscous 
timescale at ri/M = 100 is 

ivis ~ 2.54 X lO'"^ a^Hh^^Mg sec. (37) 

Hence for the case Mq — l^rriE ~ 0.1 and a ~ 0.1 we 
obtain tvis = 8.1 yr. For other cases, ri may he in the 
gas-pressure dominated region and the estimate for tvis 
will be quite different (29j . 

In reality, secular gas inflow in the hollow following 
BHBH merger is likely controlled by magnetic fields that 
are driven turbulent by MRI and other instabilities [49| . 
The magnetic field thereby acts like an effective turbulent 
viscosity, but such a viscosity might exhibit very different 
scaling behavior from the standard a— disk model. This 
issue, which is the subject of considerable interest and 
current research in the case of stationary disk accretion 
onto black holes, is clearly relevant in the case of tran- 
sient, post-merger BHBH accretion, but it is beyond the 
scope of this paper and we shall not pursue it further. 

C. Super-Eddington Accretion? 

The evolution calculated above assumes that the rate 
at which gas enters the disk is determined by conditions 
at large distance from the remnant black hole and re- 
mains constant in time during the transient. In cases 
for which the ambient gas densities are sufficiently low 
the merging BHBH may be 'gas-starved' and the ac- 
cretion rate will be safely sub-Eddington during the en- 
tire evolution. However, an interesting scenario arises if 
there is sufficient gas in the neighborhood of the merg- 
ing ('gas-rich') BHBH for the luminosity to exceed the 
critical Eddington limit Le- Suppose, for example, that 
already by the time of merger the accretion rate m(0) 



at ri in the initial hoUowed-out disk is sufficiently high 
to drive the luminosity i(0) to the critical Eddington 
value Le- Then accretion at this same rate would cause 
the luminosity to exceed L e during the subsequent tran- 
sient phase, as the total disk luminosity steadily increases 
(see Fig. [SJ . According to equation this initial ac- 
cretion rate will be larger than the critical Eddington 
rate in the final equilibrium disk by a substantial fac- 
tor: m{oo) = m(0) = {ri/risco)iTLE > ttie- Ignore for 
the moment the possibility that the geometrically thin- 
disk approximation could break down as L approaches 
Le |50|. Since the viscous timescale increases rapidly 
with r (see, e.g., equation [10]) the outer disk has insuffi- 
cient time to adjust quasistatically and reduce the inflow 
rate to maintain L{t) ai Le as gas fills the hollow. What 
happens instead is unclear and may be complicated and 
highly dynamical. One possibility is that the disk puffs 
up, the accretion becomes almost spherical as the gas fills 
the hollow, and the local radiation pressure blows excess 
matter away to maintain L w Le- Alternatively, density 
inhomogeneities resulting from the photon-bubble insta- 
bility [51, 52] might keep the disk geometrically thin and 
allow radiation to escape through porous regions of very 
low density. In the latter case the escaping luminosity 
could exceed by a considerable factor (10 — 100), ac- 
companied by radiation-driven mass loss |53| . A full un- 
derstanding of this super-Eddington scenario may require 
global, general relativistic, radiation-MHD simulations in 
3-1-1 dimensions. Such a treatment may well be worth- 
while, since this transition from quiescent sub-Eddington 
accretion to highly dynamical super-Eddington flow may 
provide another observable signature of binary BHBH 
mergers in 'gas-rich' regions. 
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